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A zero range approach is used to model resonant two-body interactions between three identical bosons. A dimen- 
2 ' sionless phase parametrizes the three-body boundary condition while the scattering length enters the Bethe-Peierls 

boundary condition. The model is solved exactly at zero energy for any value of the scattering length, positive or 
negative. From this solution, an analytical expression for the rate of three-body recombination to the universal 
shallow dimer is extracted. 
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^^ ' 1. Introduction 

("^ ' The few-body problem is a fundamental tool in the study of dilute degenerate gases such as ultracold atomic 

^D . systems [1]. It provides a basis to construct effective low energy many-body theories [2,3,4,5,6,7]. In addition, two- 

and three-body scattering properties determine the recombination and loss rates [8,9,10,11] in atomic vapors and 

are thus of practical importance. The last decade has witnessed the rapid progress in the Feshbach resonance 

technique [12] which allows to tune experimentally the two-body scattering length a. An intense experimental 

^\^ ' and theoretical activity has then developed to investigate ultracold atoms with resonant short-ranged interactions 

JH , characterized by large scattering lengths [5]. Universality [1,13] emerges in this regime, in particular at the level 

of few-body physics [14], for which details of the potential at short distances become unessential. In this context, 

universal predictions can be captured by zero range potentials where atomic interaction enters only via boundary 

conditions imposed on the wavefunctions. 

For identical bosons, Efimov [15] has predicted a universal hierarchy of shallow bound trimer states. The 
corresponding energies form a geometric spectrum, which signals a discrete scaling symmetry. At unitarity, there 
is an infinite number of such trimers with an accumulation point at the continuum threshold. The first evidence 
of Efimov physics by a resonant enhancement of three-body recombination was reported in Ref. [16] with ^''"^Cs 
atoms. This observation has been followed by many experiments on Efimov physics involving various atoms, 



Email addresses: mora@lpa.ens.fr (Christophe Mora), a.gogolin@ic.ac.uk (Alexander O. Gogolin), 
egger@thphy.uni-duesseldorf.de (Reinhold Egger). 

Preprint submitted to Elsevier Science July 29, 2010 



fermionic or bosonic [17,18,19,20,21,22,23]. Note that the Efimov effect does not require identical bosons. For 
three atoms with identical mass, operators exchanging atoms commute with the Hamiltonian. The presence of 
a fully symmetric sector is then sufficient for the Efimov effect to develop. This rules out for instance identical 
fermions but it is compatible with fermions in three different spin states. 

A vast literature is already devoted to the study of the three-body problem. There exist recent and compre- 
hensive reviews [1,13,24,25] on the subject covering a variety of approximate and exact techniques developed for 
this problem. In a remarkable series of papers, Gasaneo, Macek and Ovchinnikov have introduced the so-called 
Sturmian method which allows to derive some exact results in the universal, or zero range, case. They first solved 
a simple model [26] where only two atoms interact, and then obtained an exact solution [27] for the three-body 
problem at zero negative energy and for the atom-dimer scattering length. Finally, they were able to extend [28] 
the exact solution to an arbitrary positive energy. The zero-energy limit is feasible but remains quite involved. 

For negative energies E < and a positive scattering length, the system of three bosons only exists as a single 
atom on one side, and a dimer on the other side binding the two remaining atoms. This atom-dimer compound is 
entangled at short distance and decouples at large distance. The situation changes dramatically when the zero- 
energy E = threshold is crossed and a new channel opens: the three atoms can also separate freely at large 
distance. Therefore the zero-energy wavefunction is singular: it depends on whether the zero-energy limit is taken 
from positive or negative values. 

In this paper we investigate the case of three identical bosons interacting via zero range potential with an 
arbitrary scattering length a. In order to proceed with a well-defined model, the Bethe-Peierls two-body boundary 
condition is supplemented by a three-body boundary condition. The resulting model is solved exactly at zero 
positive energy. As an application, an exact formula is derived for the three-body recombination rate in agree- 
ment with previous works. The approach developed in this paper has similarities with the Sturmian method of 
Refs. [26,27,28]. It is nonetheless based on an integral equation that derives from the Schrodinger equation, a 
method pioneered in Ref. [29] . Part of the results presented in this paper have already been shown in a preceding 
letter [30]. The zero- range model is introduced and reviewed in Sec. 2. The solution of the three-boson problem at 
zero energy is explicited in Sec. 3 both for positive and negative scattering lengths. It is used in Sec. 4 to derive 
the rate for three-body recombination. Sec. 5 summarizes the results obtained in this paper. 



2. Zero-range model 

2.1. Integral equation 

We introduce a zero-range model for three identical bosons of mass m. The three-boson wavefunction is simply 
the solution of the free Schrodinger equation with, in addition, the Bethe-Peierls boundary condition 

which recovers the two-body physics, y is the distance between two bosons (denoted 1 and 2) and \/3x/2 the 
distance between their centroid and the third boson (denoted 3), the center of mass being decoupled. The two 
other two-body boundary conditions follow from symmetrization of the wavefunction. The Schrodinger equation 
on tl) can be transformed into an integral equation for the reduced atom-dimer function /(x). The corresponding 
procedure is now standard [6,7,11,29,31,32,33] and we will only briefly review it here. Applying the Laplacian on 
the three boundary conditions implied from Eq. (1), the Schrodinger equation acquires a source term 

V^ + V^ + ^)v(x,y)=5(x,y) (2) 

with 5(x,y) = (1 + (5)/(x)5(y). The operator Q = P23 -I- A3, where Pij exchanges atoms i and j, ensures the 
proper bosonic symmetry. 



Using the Green's function 



G£(pi,P2) = ^-- — 2 ' (3) 



the solution of Eq. (2) takes the form 

^p{x,y) = ipo{x,y)e(E) + / d^a;' d''y'G„_B/ft2(x - x',y - y')5(x',y'), (4) 



where tpo, a symmetrized combination of plane waves, describes the incoming wave, -i/'o is solution of the homoge- 
neous part of Eq. (2). 6{E) is the Heavyside function. The integral equation is closed on the function / by taking 
the J/ — >■ limit in Eq. (4) with the Bethe-Peierls boundary condition (1). This results in 

(t^E/n^ - i) /(x) = 4^^o(x, 0) e{E), (5) 

where the expression of the operator L^ is given in momentum space by 

2.2. Three-body boundary condition 

The integral equation (5) was first derived by Skorniakov and Ter-Martirosian [29], and shown later by Danilov [34] 
to be ill-defined with a dense and unbounded spectrum at negative energy. This unphysical prediction is related to 
the Thomas collapse[35], or fall to the center effect, at vanishing distances between the three bosons. It suggests 
that an additional boundary condition [36,37] is necessary if we wish to define a proper zero-range model for 
this problem. The operator L^ does not break the rotational symmetry and thus decouples the different angular 
momenta. We shall concentrate on the rotationally invariant s-wave channel, /(x) = ,f{x), where the fall to the 
center pathology does occur. At small distances a; — > 0, the inverse scattering length 1/a, the energy e and the in- 
coming wavefunction xpo can be neglected in Eq.(5). Therefore all wavefunctions become solutions of the universal 
homogeneous equation 

Uf{x) = 0. (7) 

The operator Lo has a simple action on power functions (see appendix A), 

Zo.-<-^' = 0(.).-<-^) 0M = ^^^fi^), (8) 

_ 8^inW6) (g) 

^vcos{v'k/2) 

has been defined. Giy) has an infinite number of roots on the real axis and two complex conjugate imaginary roots 

at v = ±iso with so — 1.00624. They correspond to the incoming and outgoing solutions to Eq. (7), , that 

dominate the a:: — > asymptotic behavior of /. All other solutions to Eq. (7) attenuate more rapidly as x tends to 
zero. The limiting form of f{x) is thus not fully determined since any combination of these two complex conjugate 
solutions is admissible. This has been shown [34] to lead to an unphysical dense and unbounded spectrum: there 
exists an eigenstate for any (negative) value of the energy E. 

The three-body problem becomes well-defined if one imposes the additional boundary condition [15,34,36] 

/WU^oOc ""[^°^"^;/^^ + ^' . (10) 

for the asymptotic behaviour at vanishing x, i.e. when the three bosons all coincide. The length (. does not 
bear any physical meaning, its only purpose is to set up a reference length in the logarithm. It is defined up 
to a multiplicative factor e'^'"" and can therefore take arbitrarily small or large values. All possible boundary 
conditions are exhausted when (fi is taken in the interval [0, tt], and changing the value of £ only amounts to shift 
the origin of (fi. Note that fi has a clear physical meaning: it can be thought of as the phase shift between the 
incoming and outgoing waves with respect to the origin a; = 0. The amplitudes of these two waves are equal as 
a result of flux conservation. This property is however violated once the recombination to deep dimer states is 
included in the formalism [38,39]. 

Eq. (10) can be given a transparent physical meaning. The phase shift fi is determined by short-range three-body 
physics occurring at inter-particles distances for which the zero-range approach is not valid. This short-distance 
physics is not sensitive to the larger length scales of the problem - such as the scattering length a - and fi} appears 
as a universal phase that applies to all eigenstates of the three-body problem, fi thus plays the role of a three- 
body parameter. It is the analogue of the scattering length a in the two-body problem. Finally, Eq. (10) exhibits 
a discrete scaling symmetry x — > xe'^'''° from which the geometrical pattern of Efimov states originates. 



where the function 



2.3. Two-channel model 

For a model with a finite range potential of size lo, the scaling limit is such that lo is sent to zero while 
the scattering length a is kept finite [1]. A zero range model is expected to describe directly the scaling limit 
of a more realistic finite range potential. Enforcing a particular set of boundary conditions to the solutions of 
the free Schrodinger equation does not automatically constitute a well-defined zero range model. It is indeed 
notoriously difficult to prove rigorously that the corresponding model is self-adjoint. We have circumvented these 
difficulties in Ref. [30] by introducing a two-channel model: interactions between atoms in an open channel are 
mediated by a molecular state in a closed channel. In addition to the scattering length a, this model introduces an 
additional length R* that corresponds to the size of the molecular state coupled to the open channel continuum. 
Similarly to Refs. [11,40], it amounts to give some energy dependence to the scattering length a, and the length 
K* is proportional to the s-wave effective range [1]. In addition to providing a quantitative description of a 
narrow Feshbach resonance, the two-channel model can also be seen as a toy model describing a generic two-body 
interaction in the scaling limit. It is moreover self-adjoint by construction and thus regularizes the pathologies 
associated with the three-body problem. 

The scaling limit _R* — > of the two-channel model introduced in Ref. [30] can be checked to reproduce the 
boundary conditions Eqs. (1) and (10). Within the two-channel model with K* — > 0, Eq. (5) is replaced by 

(-R-Vl + Le--\ /(x) = 4^7/,o(x, 0) e{E), (11) 

where e — mE/lf' . Since R* <C a, this equation can be split into two integral equations. For x ^ R* , the term oc R* 
inside the parenthesis can be neglected and /(x) is solution of Eq. (5). On the other hand, for a; <gC a, l/\/J£|, all 
a and e dependences disappear and Eq. (11) simplifies to 

R*\/i + io) /(x) = 0. (12) 

The complete solutions to Eq. (11) are obtained by matching the asymptotes of the different solutions in the 
region R* <^ x <^ a, where the validity domains of Eqs. (5) and (12) overlap and Eq. (7) holds. In other words, 
Eq. (12) serves as a boundary condition for the solutions of Eq. (5). As detailed in appendix A, the solution to 
Eq. (12) is obtained following Ref. [30]. For R* <^ x <^ a, one recovers the asymptotic form (10) with Ro = 
£e-^/^« ~0.577R*. 



3. Zero energy solution 

Having established a well-defined zero-range approach for three bosons, we proceed with the construction of 
the solution at zero (positive) energy. We shall solve Eq. (5) at £ = 0+ with the boundary condition (10). 
The incoming wave is generally a symmetrized combination of plane waves, 

where the summation includes all permutations P„ of {1, 2, 3}. V is the system volume and the wavefunction has 
been normalized to one. All wavevectors ki vanish at zero energy, hence i/)o(x,0) = yG/V^. Let us define the 
auxiliary function f{x) such that 

fix) = -4^a^/4 (^-^^ + ^'^(0)/») , (13) 



where 4>{u) is defined Eq. (8). In particular, one has ^(-1) = -4/\/3 and <j)(0) = |G(0) = |[1 - 47r/(3V3)]. The 
property (8) of the operator Lo on power functions yields the following integral equation 

Lo - i) /(x) = J. (14) 

An exact solution to this equation is given in the next four subsections. Subsections 3.1, 3.2 and 3.3 are devoted 
to the case of a positive scattering length a > 0, while the case a < O is also exactly solved in Subsection 3.4. An 
alternative way of expressing the solutions to Eq. (14) has been proposed in Ref. [41]. 



3.1. Solutions of the homogeneous equation 

We consider the case a > 0. The homogeneous equation associated to Eq. (14) corresponds to replace the 
right-hand-side of Eq. (14) by zero. The solution to this integral equation can be written as a Barnes-type contour 
integral 

Ptix) = - r°° ^C±{v)V(u)sM^u/2) (-Y\ (15) 

with the Gamma function V{z). The integration contour runs on the right of the imaginary axis with a small 
and positive real part. C± is a set of two functions to be determined. We apply Lo to Eq. (15), and then use the 
identities fTiy) = V{i> -\- 1), cos{-ki//2) = sin[7r(!/ + l)/2] and Eq. (8) to obtain 

LoPH^)^- — G(M)C±(i.)r(M+l)sin[^(i.+ l)/2](-) . (16) 

■^ J —ice ^*^ ^'^' 

Assuming the identity 

C±{jy + l)^G{i^)C±{v) (17) 

and the property that C±{i>) has no singularities (poles) within the stripe < Re(i^) < 2 of the complex plane, 
the integral contour can be shifted back, v -^ u — 1, with the required result Lqj3^{x) — I3^{x)/a. Using the 
Weierstrass theorem, G{i') can be expressed in terms of its poles ±6p, where bp = 2p + 1 (p G N) and its zeros 
±Up in the complex plane, namely 

+00 2 2 

Gw - n ^- (18) 

There is an infinite number of zeros on the real axis, with ui = 4, U2 — 4.6 . . ., and exactly two on the imaginary 
axis, Mo = ±iso. For p ^ 1, we find the analytical estimate 

8 [2, mod(2p-hl,6) = 3. 

Up ^bp~\ — X < (19) 

7rV3(2p-fl) [^ -1, mod(2p + l,6) = 1,5. 

Using the representation (18) and the aforementioned identity !^r(i^) = r(!/ + 1), it can be verified that the two 
functions 

C±(.) = f[ l^r^^lrn^'^V'^l (20) 

pJl) r(i^ + bp)r(i -ly + up) 

with K() = ±iso respectively, are solutions of the identity (17) with no singularity in the stripe < Re(i^) < 2. To 
be more specific, the closest singularities from the forbidden stripe are for u = 2 and v = =FJso for the functions 
C±{u) respectively. One also verifies that the integral in Eq. (15) is well-defined with Eq. (20), i.e. the integrand is 
integrable. In fact, it is physically expected that the two functions (20) exhaust the set of functions that verify the 
requirements: (i) integrability, (ii) absence of poles in the forbidden stripe and, (iii) solution of the identity (17). 
Remarkably, Ref. [27] has introduced an alternative expression for the function C±{i^). The equivalence with the 
definition used here is not straightforward to show. 

To summarize our findings, we have exhibited two independent solutions of the homogeneous equation corre- 
sponding to Eq. (14) for a > 0. They are given by Eq. (15) with Eq. (20). 

3.2. Asymptotic behaviour at sm,all x for I3f{x) 

As we shall see below, the solutions I3i{x) give rise to the incoming and outgoing asymptotes e "'"^/a:: as 
a; — > 0. The proper combination of fi^ and fi^ is therefore necessary to match the boundary condition (10). 

The asymptotic behaviour for x <^ ais obtained by closing the integration contour on the left with a semi-circle 
at infinity in the left half plane. The pole ai v — =Fiso gives the dominant contribution with the result 

P+{x) ~ ^sinh(^) |Res(C+,-iso)r(-Jso)| 6'''°+^° '"("/"') oc ^^^^, p^ (x) = [/3+(a;)]* (21) 

tx \ 2 / X 

where So is the phase of Kes{C+,—iso)T{—iso) and Res(C+,— iso) denotes the residue of C+{i^) aX v = —iso. 
Using the result (B.2) of appendix B, the phase So can be written So = n^ + Argr(— iso) — 1.5875 [27]. Here 7 is 

7 = -i - -ArgC+(iso) ~ -0.090518155, (22) 

2 TV 



correcting a typo in Eq. (17) of Ref. [30]. The linear combination /3i = e "'i/3+ + e""^^ |3-^ , with 

5i = So + So ln(i?o/ffl) = 7r7 + Argr(-iso) + So ln(i?o/a), (23) 

and Ro = te^"^'"" , recovers the three-body boundary condition (10). 

3.3. Complete solution for a > 

A solution to Eq. (14) is derived using a Barnes- type integral similar to the homogeneous case. It reads 

132 (x) = B (e-'*i /32+ (x) - e'^^ /S^ (x)) , (24) 



with 



The phase Si is defined Eq. (23), and 



B= , (25) 

2icos5i^-G(0) 



/?rW = -/ '°°|^c±Wr(.)^j^^p)"\ (26) 

^ J -too 2i7r tan(7ri^) Va/ 

with the functions C±(f) given Eq. (20). Before expliciting this solution, it can be checked that the boundary 
condition (10) is reproduced in the limit x <^ a. Similarly to Sec. 3.2, the integral contour in Eq. (26) can be closed 
around the left half plane. The asymptotic form is dominated by the closest poles located at v = ±iso and p = Q. 
The contribution from the v = pole is easily computed. Using the property C±{v)C±{l — v) = 1 and Eq. (B.5) 
derived in appendix B, it is shown to cancel exactly the term oc 1/x in Eq. (13) when f(x) is expressed in terms of 
I32{x). As a result, the asymptotic form for f{x) is solely determined by the poles At v = ±Jso. Following Sec. 3.2, 
the contributions of these two poles recover the behaviour (10). 
In order to verify Eq. (24), we apply Lo to Eq. (26). The result is 

t r.±, . 1 /'+'°°+^ du ^ , , , ,sin(7riy/2) ( x\-" , , 

xJ^ioo+i 2i7r tan(7ri/) \aJ 

where the contour lies slightly on the right of the line Re(i^) = 1. The main difi'erence with the homogeneous case 
of Sec. 3.1 is that the integrand now has a pole in the stripe < Re(!/) < 2. This pole is located at i^ = 1 due to 
the tan(7rz^) in the denominator. The integral contour can nevertheless be shifted back, v -^ v — 1, with the result 
/Jj {x)/a and an additional contribution that can be evaluated from the pole's residue at i^ = 1, namely 

Lo/?,±(x) = ^ + ^±ili4, (28) 

a TV X'' 

where C±(l) = ±i^—G{0) from Eq. (B.5) in appendix B. Using Eq. (28), the function I32{x) from Eq. (24) is 
found to be solution of the integral Eq. (14). 

To summarize, the general solution to the integral equation (14) is of the form f{x) = Ii2{x) +C/3i(a:), where C 
is an arbitrary complex coefficient. 

3.4. Exact solution for o < 

We finally discuss the case of negative scattering length a < 0. In that case, the operator Lo + l/l^l does not 
have a homogeneous solution and can be inverted. The physical meaning of this property is that (universal) dimers 
do not form for a < 0. The atom-dimer sector is therefore absent at large distances between the atoms. Eq. (14), 
with the boundary condition (10), possesses a unique solution, namely 

P^ix) = B {e-^''Pt{x) - e^'^p^ix)) , (29) 

with B given Eq. (25) and 



**w-^/,:£-*<^w^>=£?(wi ^ *») 



Similarly to Sec. (3.3), the boundary condition (10) is verified by closing the contour around the left hand plane 
and evaluating the residues at i/ = and v — ±iso. Applying Lq to /Js amounts to shift the integral contour 
v ^ V + 1 with a sign change. A pole at i^ = 1 is encountered when the contour is shifted back, leading to 

Lol3,{x)--^-——-^. (31) 

This result, inserted into Eq. (29), shows that P^ix) is solution of Eq. (14) for a < 0. 



4. Rate of three-body recombination 

Three-body recombination is a collision process in which three incoming atoms form a (dimer) molecule and an 
atom. The binding energy of the dimer is converted into kinetic energy, and recombination thus produces losses 
of atoms from the trap in ultracold atomic vapors. This effect is often important as it can limit experimentally 
the lifetime of quantum gases. It also offers a convenient experimental tool to probe resonances in the (few)three- 
body problem that are usually accompanied by increasing recombination [42] and therefore loss rate. In the 
case of resonant interactions (large a > 0), the shallow dimer has universal features and recombination can be 
quantitatively captured by a zero-range model such as the one discussed in this paper. 

4.1. Formalism and result for the recombination rate 

The general solution f{x) — l32{x) -f C/3i(a;), that was derived in Sec. 3.3, can be inserted in Eq. (13) to find 
f{x), which gives access to the complete three-boson wavefunction -i/;(x,y) via Eq. (4). The resulting expression 
is cumbersome and we shall not write it here. Instead, the asymptotic properties of tp can be discussed quite 
generally. Two relevant sectors emerge upon considering the large distance asymptotes. Sector I: for x — > +oo and 
y — > -l-oo, the wavefunction corresponds to three free atoms in a symmetrized combination of plane waves. Sector 
II: the second sector splits itself into three equivalent domains. Taking x — > +oo with y fixed describes the atom 
3 and a dimer, formed by atoms 1 and 2, flying apart. The two remaining domains are obtained by exchanging 
the single atom 3 by 1 or 2. 

For both sectors, there are incoming and outgoing waves whose relative coefficients depend on the constant 
C. In order to calculate the three-body recombination rate, the following scattering situation is considered: the 
incoming wave is formed solely by three free atoms and the flux probability to leak into the three domains of the 
atom-dimer sector is computed. The constant C is therefore chosen in order to cancel the atom-dimer incoming 
wave. 

We use the notation r = y/Sx/2 for the distance between the atom 3 and the 1-2 centroid. For x — > +oo with y 
fixed, the asymptotical form ip{x,y) ~ <^o(y)0ad(f) is expected, where cpoiy) ~ e^^''^/{y\/2-Ka) is the normalized 
two-body bound state (dimer) wavefunction, and 

^.4r)=A^--^—, (32) 

where ^ is a constant that will be extracted from the exact solution. The wavefunction 0ad(?') describes the atom- 
dimer relative motion corresponding to an outgoing scattered wave by the recombination process. The J/ — >■ limit 
in the asymptotical form of V'(x,y) can be compared to the Bethe-Peierls boundary condition (1) to extract / 
with the result 

ix/a 

f{x)^A- , (33) 

X 

in agreement with the behaviour of f{x) at large x that shall be derived in Sec. 4.2. Interestingly, this result shows 
that the large x study of f{x) is sufficient to determine the wavefunction in the atom-dimer sector through the 
knowledge of A. A similar reasoning also indicates that a vanishing atom-dimer incoming wave is equivalent to a 
vanishing e~^^'°' jx term in j{x) at large x. 

Eq. (32) can be used to compute the rate of three-body recombination. The current probability associated to 
(?!'ad(r) is given by j[r) = 2i(2m/3) [</'ad(^)9r</'ad(j") — 9r0ad('")0ad(»")] wherc 2m/3 is the atom-dimer reduced mass. 
The probability flux "l>i to leave the scattering region in one of the atom-dimer domains (i denotes the indices of 



the single atom) can be obtained by integrating j{r) over a sphere with an infinitely large radius. An additional 
integration over the free center of mass position multiplies by the volume V and one finds 

- = $, = V hm 47rO r = -^ A^, 34 

where $ denotes the total flux to the atom-dimer sector, or sector II. It is also possible to relate this flux to the 
large x form of the auxiliary function / deflned Eq. (13), namely 

ix/a 

f{x)^A , (35a) 

X 

<&, = 192^/37^'aV'(0)-i^, (35b) 

m V ^ 

where (f){p) is defined Eq. (8) and 0(0) = |G(0) = |[1 — 47r/(3\/3)]. The value of A can be extracted from the 
exact solution of Eq. (14) derived in Sec. 3.3. This calculation is postponed to Sec. 4.2. 

For a dilute gas of A'^ bosons with density n — N/V, the three-body recombination rate Qrcc is defined [1] such 
that the number of recombination events per time and per volume is aiccn'^- Hence, 

[NV 1 diV,ec-cvcnts <S> N {N - 1){N - 2) ^N'' .^^, 

""-[vj =V di =V 6 ""V— (^^^ 

where — - — ^-2^ — ^-^ is the number of triplets among the N bosons. The resulting expression for the three-body 
recombination rate is 

Qrcc = 96^/37r^aV^(0) — 1^1^ (37) 

4.2. Wavefunction at large x 

For X ^ a, the integrands in the expressions of /3i(a;) and P2{x), Eqs. (15) and (24), develop rapid oscillations. 
The asymptotical behaviour can thus be deduced from a saddle-point analysis that is detailed in appendix C. 
Focusing on the exact solution, 

f{x) = P2{x)+ C/3i {x) = B (e-'*i pt (x) - e^'^ p^ {x)) + C (e"'*! /?+ {x) + e^'^ ft" (x)) , (38) 

the incoming atom-dimer wave (oc e^'^'"/a;) is seen to vanish with the choice 

C=- ' \ \.^ . (39) 

Inserting this result back into the large x^ a expression of f{x), Eq. (35a) is recovered with the coefficient 

.C 1 _ e^^'^'oe^"'! 
^B l + e-2''«oe2»'5i 

Some simple algebraic manipulations further lead to the expression 



_4 ^ ^^^iS^^S,)ae^ 1^ _ ^-2.s,^2^8, 



1 



(40) 



„KS^-Si)^a sin(5i cos Si 



A^Be"-'^-"^' — ' \; ■ (41) 



The identity 

(42) 



1 - e-2^'i 



1 1 g — 27rsog — 2i(5i 1 g — 27rso * 

with the phase 5r = di — Arg(l + e^'^'^""e^' ^) is used together with the value of B from Eq. (25) to derive 

r IT a jtg _A 1 sincJr /,„N 

^-G(0) smh(7rso) 

The rate for three-body recombination is finally given by 

_ 12871^ (47r - 3^3) fi^a'^ 



"rcc — -TT-, ^ sin (Sr). (44) 

smh (ttso) rn 

This function shows periodic oscillations as a function of ln(i?o/a), with a period given by the scaling factor e^"" , 
and a dimensionless amplitude 1287r^ (47r — 3\/3)/sinh^(7rso) — 67.1177. In practice, e~'^^"° ~ 0.0012 is a very 
small number such that Si ~ Sr and the oscillations are almost sinusoidal. 



Instead of Eq. (44), an alternative formula can be obtained for the three-body recombination rate by taking 
directly the absolute value squared of Eq. (41) and using that le"^" + g-'^sog-ai^ip ^ 4sinh^(7rso) + 4cos^(<5i). 
The result reads 

_ 12871^ (47r - 3^/3) ft^a^ sin2(5i) sinh^(7rso) 

sinh^(7rso) m sinh^(7rso) + cos2(5i) ' 

which coincides with Eq. (44) and recovers the result of Refs. [28,30,43,41]. Ref. [13] has pointed out that the 
expression (44) also derives from the use of the optical theorem and the knowledge of the elastic boson-diboson 
scattering solution at zero (negative) energy determined in Ref. [27]. 



5. Conclusion 

We have solved exactly the wavefunction at vanishing energy of three identical bosons with zero range interac- 
tions. The zero range model requires a set of two boundary conditions. The first one is the standard Bethe-Peierls 
condition when two atoms meet. An additional three-body boundary condition when the three atoms coalesce is 
then necessary. It is due to the fall to the center effect, also responsible for the emergence of the Efimov bound 
states. It is characterized by a dimensionless parameter (p which describes the phase shift between the outgoing 
and the incoming waves toward the region where the three atoms coincide. Close to this region, the wavefunc- 
tion displays oscillation as a function of ln{R/a) where R denotes the hyperradius of the three bosons. These 
logarithmic oscillations thus appear in most observables and in particular in the three-body recombination rate. 

The derivation of the exact solution is facilitated by the fact that the operator L^ does not involve any energy 
scale at zero energy: it is homogeneous. As a result, its action is stable in the space of power functions. The trick 
is then to write the exact solution as a Barnes-type contour integral over a power function where the power is 
the integration variable. Acting with Lq simply amounts to shift the integration contour. The integrand has to be 
chosen such that no pole is crossed when the contour is shifted back to its original position. 

The solution is unique in the case of a negative scattering length where no dimer can be formed. This is in 
contrast with the case of a positive scattering length where the solution is not unique but can be parametrized 
by a single parameter. It depends on the balance between the atom-dimer and the three-particle sectors in the 
incoming wave. The three-body recombination rate is evaluated by canceling the incoming atom-dimer wave and 
then by calculating the prefactor of the outgoing atom-dimer wave. The result is an exact analytical expression 
that exhibits the expected logarithmic oscillations. 



Appendix A. Boundary condition from the two-channel model 

It is more convenient to solve Eq. (12) in Fourier space, or 

(R*k'+Lo)f{]i)^0. (A.l) 



The corresponding operator is rotationally invariant and thus decouples the different partial waves. For all angular 
momenta I > 0, the limit i?* — > can be taken directly in Eq. (A.l), and coincides with the zero range approach. 
The s-wave sector I = 0, with /(k) = f{k), requires a particular treatment in the scaling limit R* — >■ 0. One 
notices that the operator Lq is homogeneous in the sense that it does not involve any scale. It implies that Lq 
acts simply on power functions, namely 

Lofc"""' = G(//)fc-"-\ (A.2) 

where G{u) is given by Eq. (9). This result is obtained with the change of variable k' — ke^ in Eq. (6) and the 
integral 

+°° dC ^i-e /^ e^^ + e^ + 1 ^ _ 1 sm{uTv/6) 



2tt^ Ve2«-e«-hiy ucos{vk/2)' ^ ' -* 

Using the property Eq. (A.2), one verifies that the Barnes- type contour integral (on the right of the imaginary 
axis) 

— C{u){kR')-'-\ (A.4) 



is solution of Eq. (A.l). The function 

^ ' sm[7r(i/ - iSo)] sin[7r(i/ + iso)] ^ ' ' 

such that G{iy)C{i') — — C(m + 1), has no poles within the stripe < Re(i/) < 2 of the complex plane. Applying 
L() on Eq. (A. 4) and shifting the contour as i^ — >■ i/ — 1 indeed recovers Eq. (A.l). 

The large x ^ R* behaviour of f{x) is obtained from the asymptote of f{k) at small kR* <^ 1. By closing the 
integration contour around the left half plane in Eq. (A. 4), f{k) is evaluated as a sum over terms with powers at 
least higher than k^^. The poles aX u = ±iso give the dominant contribution at small k. After Fourier transform, 
the result reads 

^^^,^^sin[.oln(^/i^o)l^ (A.6) 

with Ro = R' exp [(7r(7 + 1/2) + Argr(iso))/so] ~ 0.577i?*, where 7 is given Eq. (22). 



Appendix B. Useful formulas 

We here list some useful relations. 

First, the residue of C+ at i/ = —iso is given by 

Res(C+, ~i.o) = , r(2 + i.o) ^ ^ n-^s, + u,)nl+is,+b,) 

where the p = terms have been singled out from the convergent infinite product. The general identity r(2)r(l — 
z) = 7r/sin(7rz) is used at 2; = — 2iso to obtain 

Res(C+, -is„) = '^"^^^'"''^ [C+(iso)]*, (B.2) 

ITT 

with the help of the property r(2:*) = [r(z)]*. 

Second, using the definition (20) of the functions C±, we can write 

On the other hand, the product representation Eq. (18) of G{v) can be evaluated at i^ = leading to 



with 1 - 47r/(3^/3) ~ -1.41839. As a result, the identity 



C±(1)=±V-G(0), (B.5) 

is obtained. 



Appendix C. Saddle-point analysis at large x 



The integral expressions (15) and (24), for /3j and fi^, can be evaluated at large a; 3> a. In contrast with the 
small X case, the integration contour can not be closed on the right half plane. Nevertheless, the integrand exhibits 
rapid oscillations along the imaginary axis as x is increased. A saddle-point, or stationary phase, approximation 
is thus carried out to extract the asymptotical behaviour. 

For large i^ 3> 1, we have \i\V(v) ~ In yphi + {v — 1/2) \av — v, and the coefficient C^{v) defined in Eq. (20) 

becomes 

+00 

lnC+(j.) ~ {^{v) - v|/(l - v)\ Y^iu, - b,), (C.l) 
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with the digamma function ^(i^) — r'(i')/r(i^). Equation (19) can be used to check that the summation over p 
indeed converges. Since ^(i^) ~ Ini^ for u ^ 1, the hmits 

hm C+(!/) = eT(*^°°+""«^ (C.2) 

are obtained with the phase S^o = ^''^^t=ii'^p ^ ^p) + "" ~ 1.736 in agreement with Rcf. [27]. Noting that 
sin(7ri^/2) ~ ±ie^*''''"/2 for u — > ±ioo, we rescale the integral with u = xz and approximate 

^^~eA^s^+.so)^.a(z)^ (C.3) 

-ioo VZ 

We have defined Sz = sgn (Imz) and the function 



9i^) 



In 2 — 1 + In a - 



(C.4) 



The integral (C.3) is dominated by two symmetric saddle-points on the imaginary axis, Zq — ±i/a, derived from 
the condition g'{zo) = 0. Moreover, one has gizg) = ^i/a and (/"(zg ) = ^ia. After integration over the quadratic 
fiuctuations around the two saddle-points, the asymptotic form 

/3t{x) ~ -i^ ( e"'«e"^~ ^— - £-''»« e-'*~ — — j (C.5) 

is obtained for x ^ a. The same result is derived for 13^ (x) with so replaced by —so. In this asymptotic region 
X ^ a, Eq. (C.5) describes a combination of incoming and outgoing atom-dimer waves. 

The calculation is similar for the case of /^j • The difference is expressed by the following limit 



sin(7ri//2) 
tan(7ri^) 

for I' — > ±ioo. The final result reads 



(C.6) 



P+{x) ~ ^ e"'V°° + e-"'«e-'''~ , (C.7) 

2 \ X ^ J 

with So — > —So for /?^(a::). Again incoming and outgoing waves coexist in this asymptotic form. 
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